Can tetraneutron exist from theoretical point of view? 
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Abstract 

A theoretical possibility is shown for the bound state of a tetraneutron to exist in the case of 
the proposed neutron-neutron potentials in the singlet state with two attractive wells separated 
by a repulsive barrier. The anomalous behaviours are revealed for the calculated size, density 
distribution, and pair correlation functions of a hypothetical tetraneutron. 
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1 Introduction 



^ \ The mysterious fact of the experimental registration of tetraneutrons P^E] in a reaction with loosely 
^ ' bound radioactive ^^Be renewed the attention to theoretical attempts of understanding the problem 
^ ■ of hypothetical bound neutron systems. The experiment contradicts the rather old estimates showing 
^ . the impossibility to form a bound state of a few neutrons interacting by the standard nuclear forces 
(see review Recent attempts to study this problem using more modern methods of calculation 
of four-particle systems |U El E! also indicate the impossibility of binding the four-neutron system 
without adding some exotic many-particle interaction potentials. Moreover, the absence of resonances 
in a four- neutron system with standard potentials was shown in [Zj. The experimental search for 
resonances in and systems |S] using other nuclear reactions had also no success. 

In the present paper, we study the theoretical problem of the possible existence of a tetraneutron 
by developing the precise methods of calculation of loosely bound states of four Fermi particles. We 
propose a special idea of constructing the neutron-neutron potentials allowing the bound tetraneutron 
to exist and simultaneously describing the standard low-energy neutron-neutron data. 



2 Basic equations 

To study the properties of the four-neutron system in the state with zero spin {S = 0) and orbital 
moment (L = 0) under assumption of the central pairwise neutron-neutron interaction potentials. 
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we solve the following Schrodinger equation for one spatial component of the wave function: 



i=l i>j=l foV(14),(23) 

(ii)^(12),(34) 

The total antisymmetric wave function of the four-neutron system is expressed in terms of the 
corresponding spin and spatial components as 

vl/'^(l, 2, 3, 4) = ($'(ri, r2, rg, v,)^ - $"(ri, r2, rg, r4)e') , (2) 

where the antisymmetric $' and symmetric $" (with respect to permutations (1 ^ 2) and (3 ^ 4)) 
spatial components are 



<l>'(ri, ra, rg, = $(ri, rs, rg, r4), (3) 
<l>"(ri, r2, rg, = — (2<l>(ri, ra, rg, r4) - $(ri, rs, rg, r4)) , 

which corresponds to the Young scheme [2,2]. In Eq. P23 is the permutation operator of 

spatial coordinates, V^{rij) and Vfijij) are, respectively, the singlet interaction potential in even 
states and the triplet one in odd states. The bound states of four neutrons are studied by solving 
the Schrodinger equation for various nuclear potentials taken in the form of a superposition of 
Gaussian functions, by using the well-known variational method with the translationally invariant 
Gaussian basis antisymmetrized with respect to the permutations of particles (1 2) and (3 4), 

<l>(ri, ra, rg, v^) = exp | - ^ ujrj j , (4) 

k=l \ i>j=l J 

where N is the basis dimension, and r^j = |rj — rj|. This basis and the special schemes necessary 
to optimize the nonlinear variational parameters u'^j enable us to carry on the calculations of loosely 
bound states with desired high accuracy. 



3 Spinless interaction model 

First, consider the simplest case of potentials independent of spin ( i.e. V^~{rij) = V^{rij)) where the 
overestimated attraction in the triplet state may only promote the binding of the four-neutron system. 
We study the bound state appearance conditions varying the coupling constant of potentials of 
different forms (further, we use the dimensionless units: where rn is 

the radius of interaction, and g is the coupling constant; f{} jm = ^XA^'l'b MeV-fm^ for neutrons). 
For the potential with one Gaussian function, IJir) = — 5'exp(— r^), a bound tetraneutron exists 
below the decay threshold (4 — > 2 + 2) only for g > 5'cr(4) = 3.911, which is 1.46 times greater than 
the critical two-particle coupling constant gcrC^) = 2.684. We notice that the reliable calculations 
need basis (ji} to be about 150 functions with the optimization of nonlinear parameters. Note that 
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a trineutron, for the same potential, can be bound below the decay threshold (3^2 + 1) only 
for g > 3. 3(7^(2). From the qualitative point of view, the similar conditions of the tetraneutron 
bound state appearance take place for other traditional purely attractive potentials. Thus, under 
the condition that a dineutron is unbound, there is no possibility, because of the Pauli principle, 
for a tetraneutron (nothing to say of a trineutron) to form a bound state with traditional attractive 
potentials. An analogous conclusion can be drawn in the case of common potentials with repulsion 
at short distances between neutrons. For example, for the widely used Volkov potential, one has 
k = gcr{4:) / Qcri'^) = 1.44. Moreover, the attempt to find a better ratio k by varying the parameters 
of the two-component potentials with attraction and short-range repulsion led us only to a potential 
U{r) = 5f (1.5 exp(— (r/0.9)^) — exp(— r^)) giving rise to k about 1.27. We can assume that it is 
impossible to form a bound system also for other standard interaction potentials with attraction 
and short-range repulsion if is unbound, which is in agreement with the recent calculations [HEIISI- 
In principle, a possibihty for a bound tetraneutron to exist, under the condition of an unbound 
^n, can be realized with some exotic pairwise interaction potentials having two regions of attraction 
separated by a repulsive barrier. An external attractive potential well of greater radius is necessary, 
first of all, to fit the experimental low-energy two-neutron scattering parameters in the singlet state, 
and it has to be in the typical range of nuclear forces of about or greater 1.5 fm. Fitting the singlet 
interaction potential, we use the following low-energy neutron-neutron scattering parameters: the 
scattering length as(nn) = —18.9 fm and effective radius ros(nn) = 2.75 fm. The internal attractive 
potential well of smaller radius is important for binding the system, while the repulsive barrier 
between the attractive wells makes the two regimes of attraction somewhat independent. In a tetra- 
neutron, the number of pairs of particles in the singlet state, as well as that in the triplet one, equals 
three, the Pauli principle reveals itself only in the triplet state, and the internal potential well acting 
in the singlet state plays the main role in binding the '^n system. A class of potentials with two 
attractive wells of different radii, which give rise to the bound state of under the assumption 
Vt~{rij) = V^{rij), is rather wide. We have a number of potentials in the form of a superposition of 
three or four Gaussian functions. One of the optimal variant is the four-component neutron-neutron 
singlet potential 

f/+(r) = g {0.43 exp (- (r/0.6)^) — exp (— r 2) + 1.085 exp (- (r/1.3)^) - 0.42 exp (- (r/1.5)^) } , 

(5) 

where the distance is measured in units of vq = 0.488519 fm. At g = gexper = 322.40, potential 
dni) reproduces the experimental values of as(nn), rQs{nn), and the commonly used recommended sin- 
glet neutron-neutron phase shift up to the energies Ei^h ~ 80 MeV. Note that there are no direct 
measurements of the neutron-neutron phase shifts. 

Fig. ^ shows the dependence of the tetraneutron energy on the coupling constant g of potential 
(0) in the "spinless" case Uf{rij) = U^{rij). Note that the decay threshold of into 2-1-2 (as well 
as the decay threshold of into 2-1-1) as a function of the coupling constant has two regimes of 
behaviour. The first regime of a rather weak binding of at g — *• gcr{'^) takes place due to the 
presence of the attraction of greater radius in potential (0) , and the second one with the almost linear 
dependence of the threshold in a wide range of coupling constants is present due to the attraction 
of smaller radius. A repulsive barrier between the attractive wells contributes to the sharpness of 
changing the two regimes of the threshold behaviour. Note also that the excited two-particle S'-state 
lies anomalously close to the ground state at g*^{2) / gcri^) = 1.12, which is caused to a great extent 
by the presence of two almost independent attractive wells in potential It is essentially important 
that, in variational calculations, the system is bound already ai g > (7^(4) = 315.2 = 0.954(yfcr(2), 
where is still unbound {g < gcri!^) = 330.42). Moreover, at the coupling constant g = gexper = 
322.40 = 0.976(7^(2), where potential (jSj) reproduces the experimental low-energy neutron-neutron 
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parameters, the system is already bound. At the same time, a trineutron with the considered 
potential is not allowed to be bound since the ratio (7cr(3)/5'cr(2) = 1.008 is greater than 1, although 
being close to it. Notice the fact that the energy dependence on g looks like almost a straight line 
parallel to the energy threshold (2 + 2) dependence in a wide interval of coupling constants {g > 1.1), 
and this line is rather close to the (2 + 2) threshold. This fact indicates that, in this region of g, a 
tetraneutron exists due to the presence of the internal potential well of smaller radius, and the 
state is of the two-dineutron cluster nature. A similar consideration concerns as well: in a wide 
interval of coupling constants, it is the cluster state (2 + 1) with the essential role of the attractive well 
of smaller radius. This is confirmed also by the approximate relation i?4„ — 2i?2„ « 2 (i?3„ — -E'2„). 
By the way, we constructed some other variants of potentials U^{r), for example, 

f/+(r) =^ {0.315 exp(-(r/0.5)^) - exp (-r^) + 1.278 exp (- (r/1.31)^) - 0.54exp (- (r/1.5)^) } , 

which yields gcr{4:) / gcrC^) = 0.9525 for to exist and reproduces the low-energy parameters of n — n 
scattering at gexper/ Qcri'^) = 0.9935 ((7^(2) = 187.5, and interaction radius tq = 1.2608 fm). Even a 
trineutron could exist for this potential upon the unbound "^n due to the ratio gcri'i) / Qcri'^) = 0.9564 
being less than unity, but the more strict condition ^^^(3) < gexper is not valid. In addition, the 
n — n singlet phase shift for this potential becomes too large already at rather low energies. We have 
found no variant of the potential obeying the condition 5'cr(3) < gexper simultaneously with giving a 
reasonable phase shift at low energies for ^n to exist in a bound state. 

It should be noted that it is necessary to carry on variational calculations with special schemes 
of optimization of basis (0)) in order to obtain the above results for '^n with reliable accuracy. In 
particular, we used about 220 functions with the optimization of the basis for this purpose. This is 
caused by both the complicated antisymmetrized four-neutron wave function of the near-threshold 
state and the potential containing essentially different components. 

4 Realistic case 

Now consider more realistic models of the triplet interaction potential Uf{r), when the conditions 
for the existence of the bound state of a tetraneutron are somewhat less appropriate. If we put 
Uf{r) to be zero in Eq. (0), we get an unbound tetraneutron within the proposed class of singlet 
potentials, because only a half of 6 pairs interacts in this case. Thus, it is necessary to have some 
additional attraction in odd orbital states to bind ^n. On the other hand, the commonly recommended 
phase shifts of the scattering in odd orbital states are rather negative corresponding to the effective 
repulsion. It appears that there exists a class of triplet potentials which together with the singlet 
potential (5) can bind the '^n system and are repulsive with the exception of the typical nuclear 
distances of about 1.5—2 fm, where they reveal some attraction correlated with the external attractive 
potential well of the singlet potential. Such a potential (in the same dimensionless units, as potential 
(jS))) may have the form 

Uf{r) = gt {2.212 exp (- (r/2)^) - 2.334 exp (- (r/sf) + exp (- (r/4)^) } (6) 

with gt = 14. This potential has negative phase shift in the P-state, although with non-monotone 
dependence. Potential (jH)) together with the singlet one (0) result in the bound state of a tetraneutron 
with the binding energy B{^n) > 0.5 MeV (this value is the variational estimation with the use of 
about 400 Gaussian functions). A more accurate calculation needs much greater efforts mainly 
because of the complicated structure of potentials and the many-component antisymmetrized wave 
function of four particles. In addition, the ultimate result for the binding energy is a few orders 
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Figure 1: Tetraneutron energy dependence on the coupling constant of potential ^ acting both in 
the singlet and triplet states (ro = 0.488519 fm is the radius of the interaction potential ©)• 
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Figure 2: Profile of tlie density distribution n(r) of neutrons of the hypothetical tetraneutron {tq is 
the same as in Fig. Q). 

of magnitude lesser than the contributions of the kinetic or potential energies calculated separately 
with rather good accuracy, and these contributions almost cancel each other having opposite signs. 
The proposed phenomenological potentials are constructed only to demonstrate the possibility for a 
tetraneutron to exist in a bound state. Moreover, one can easily change the binding energy of "^n in 
a wide range (from zero to dozens of MeV) by changing slightly potential (jHI) or potential (jSj). There 
are some reasons to assume that if these potentials with repulsive barriers are changed in such a way 
that they should not allow a tetraneutron to be bound, they may result in resonances in the system 
of four neutrons. 

Consider the main structure functions of the hypothetical tetraneutron. Note that the structure 
functions can be calculated much more accurately using basis (4) of a lesser dimension than that 
used in the calculation of the energy. Fig. El presents the one-particle density distribution of '^n 
(normalized as / n{r)dr = 1 ) versus the dimensionless distance, for potentials (0), Due to 
the Pauli principle, the density distribution has essential minimum at short distances, i.e. the 
tetraneutron is a "bubble" system with the almost Gaussian near-surface distribution of neutrons. 
The tetraneutron has anomalously small (in nuclear scale) r.m.s., (r^)^^^ = 1.704ro = 0.83 fm, which 
is caused by the attraction well of smaller radius in potential Changing the singlet potential 
one can increase the above value mainly due to an increase of the radius tq. But, in any case, the 
size of a tetraneutron, in spite of its extremely small binding energy, will be less or about the size of 
an a-particle. 

Fig- El depicts the singlet g2sij) and triplet g2tij') pair correlation functions, which reflect, to a 
great extent, the behaviour of the corresponding potentials. The singlet correlation function g2s{r) 
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Figure 3: Profiles of the singlet g2s{T), triplet g2t{r), and total g2{r) = | {g2s{r) + g2t{r)) pair corre- 
lation functions of the hypothetical bound system [tq is the same as in Fig. Q). 



has significant maximum in the region of internal short-range attraction of the singlet potential, 
since the Pauli principle does not reveal itself in the singlet state. Some decrease of g2s{f) at very 
short distances is caused by the presence of short-range repulsion in potential (0), and it should 
not be present if the repulsion were absent. The secondary maximum is present in g2s{r) due to 
the existence of the external attractive potential well of V^{r). In the triplet state, the repulsion 
at short distances makes a small contribution into the energy because of the Pauli principle, and 
the maximum of g2t{f) is located in the attractive area of the triplet potential. That is why, the 
contribution of the triplet potential to the energy of is negative, and a tetraneutron could not 
be bound without the contribution of this comparatively small effective attraction. The short-range 
attraction in the singlet state plays the main role in binding the '^n system. This is confirmed by 
calculations of the average singlet and triplet potential energy contributions. 



V^{T)g2s{r)dr+ / {r)g2t{T)dv ^ (V+> + {Vf) 



{V) 



where {V^) = —1296.7 MeV and (V^^) = —158.7 MeV, which together with the kinetic energy 

(^K^ = 1454.9 MeV result in the negative energy of the system £^4„ < —0.5 MeV indicated above. 

The total correlation function (72 (r) = | {g2s{r) + g2t{j)) reflects the average neutron pair correla- 
tions and has the main maximum at short distances and the secondary one in the region of attraction 
in the triplet state (where the singlet potential also has an external attractive well). 
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5 Conclusions 

To summarize, we note the following. 1) A tetraneutron can exist in the bound state if one assumes 
that the interaction potential in the n—n singlet state has two attractive wells separated by a repulsive 
barrier. Unfortunately, as a result, we get an anomalously high maximum in the singlet scattering 
phase shift 5s ~ 160" at the energies of neutrons of the order of 100 — 150 MeV. The problem 
of constructing the potential, which binds and gives no anomalous maximum mentioned above, 
is to be further studied. Maybe, a combination of pairwise potentials with some small intercluster 
ones can improve the situation. The nature of the assumed exotic short-range attraction is to be 
discussed as well. 2) Strange though it may seem, few-nucleon systems are essentially underbound 
with the n — n potentials (jHI) used together with the standard n—p ones. In particular, with the 
Minnesota potential used as the n — p interaction, one has the binding energies of about 6.2 MeV 
for ^H, and about 23.4 MeV for ^He. The calculated few-nucleon values may be in agreement with 
the experimental data under the condition that the fitting of the potentials is carried out taking 
into account the concordance of the regimes of attraction of all the potentials. 3) The hypothetical 
tetraneutron has abnormal small size and, at the same time, small binding energy, and this may 
serve as an additional criterion for the identification of such systems. It is interesting to study 
the probability of the tetraneutron "presence" in ^"^Be nuclei. 4) Potentials © and © satisfy the 
saturation conditions necessary for the stability of the neutron matter. Interesting and nontrivial 
questions arise concerning heavier multineutron systems with such potentials. In particular, a system 
of 8 neutrons can be a more stable system than both because of the magic number of particles, 
and because the former may be like two tetraneutron clusters. 
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